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1. Introduction 

The main aim of our paper is to develop a many-body theory approach to the calculation 
of gamma spectra from positron-atom annihilation. 

Positron annihilation has been an important tool for studying electronic and atomic 
structure of solids for over half-a-century. Its history can be traced from the first review 
of the theory of positron annihilation in solids by Ferrel (1956) to a more recent work 
by Puska and Nieminen (1994). Positron lifetimes and spectra of annihilation gamma 
quanta contain information about the structure and composition of bulk materials and 
surfaces, presence and concentration of defects and voids or pores, and their sizes, and 
electron momentum distribution. Interpretation of positron annihilation data requires 
good theoretical understanding of the process. Many-body theory was the first method 
used to determine positron annihilation rates in metals (Kahana 1963, Carbotte 1967). 
In particular, it was successful in explaining large enhancement factors that increase the 
annihilation rates above that obtained in the noninteracting electron gas approximation 
(see also Arponen and Pajanne 1979). 

In the 1990's measurements of gamma-ray spectra from positron annihilation on 
atoms and molecules in the gas phase became possible (Tang et al 1992, Coleman 
et al 1994). For He the experiment showed excellent agreement with the calculated 
spectrum obtained from an elaborate variational positron-He wavefunction (Van Rceth 
et al 1996). This work also revealed sizeable deviations of the shape spectrum from a 
Gaussian, which is often used as an approximation. For noble gas atoms of Ar, Kr and 
Xe, a careful study of the shape of the 511-keV gamma-ray hne provided an estimate of 
the contribution of positron annihilation with inner-shell electrons (Iwata et al 1997a). 
A large systematic study was conducted for a variety of inorganic molecules, alkanes, 
alkenes, aromatics and perfluorinatcd and partially fluorinated hydrocarbons (Iwata et 
al 1997b). In particular, this work determined the relative probability of annihilation 
on the fluorine atoms and on the C-H bonds for partially fluorinated hydrocarbons. 
Such information is important for achieving better understanding of very large positron 
annihilation rates on polyatomic molecules (see, e.g., Surko et al 2005). 

On the theory side, much progress in understanding the interaction of low-energy 
positrons with many-electron atoms has been obtained by using many-body theory. 
Its application to positron-atom collisions was pioneered by Amusia et al (1976) who 
considered positron scattering from He. Subsequently, the effect of virtual positronium 
(Ps) formation on positron-atom interaction was investigated for heavier noble gas 
atoms (Dzuba et al 1993, 1996). These works showed that virtual Ps formation gives 
a large contribution to the positron-atom attraction|. It leads to the emergence of 
positron-atom virtual and bound states (Dzuba et al 1995, 1996), and completely 

I This explained the success of earher polarized-orbital calculations in describing positron elastic 
scattering and annihilation for noble gas atoms (McEachran et al 1980 and references therein). In 
that approximation the positron was treated as a heavy particle, and the strength of the positron-atom 
polarization potential was overestimated, making up for the complete neglect of virtual Ps formation. 
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alters the picture of low-energy positron-atom scattering (Gribakin and King 1996). 
These papers employed an approximate method of accounting for virtual Ps formation. 
Recently, a consistent many-body theory approach has been developed. It is based on 
the summation of the electron-positron ladder diagram series (Ludlow 2003, Gribakin 
and Ludlow 2004). Its contribution is especially prominent in the calculation of 
the annihilation rate. These and earlier calculations (Dzuba et al 1993, 1996) have 
demonstrated that electron-positron correlation effects can enhance the positron-atom 
annihilation rate as much as 10^ times. 

In the paper by Iwata et al (1997a), the fraction of annihilation with inner-shell 
electrons was derived by fitting the experimental data with a linear combination of the 
gamma spectra for the valence and inner shells. These spectra were calculated in the 
simplest approximation, by using the positron wavefunction in the static field of the 
atom. In spite of a complete neglect of electron-positron correlations, the shapes of 
the measured gamma spectra were described reasonably well. On the other hand, some 
discrepancies between the theoretical and experimental spectra were obvious, adding 
uncertainty to the estimates of the inner-shell annihilation fractions. For example, the 
data for Ar was compatible with a zero contribution of annihilation with 2s and 2p 
electrons. 

Since correlations effects play such a large role in positron-atom interactions, they 
must be included in the proper theory of annihilation gamma spectra. In what follows we 
briefly recount the main facts and formulae concerning positron annihilation and gamma 
spectra using the formalism of creation-destruction operators (section Ej). We then 
proceed to derive the many-body diagrammatic expansion of the annihilation amplitude, 
and consider the significance of various terms and their relation to the total annihilation 
rate (sectional). In section |3J expressions for the 0th and Ist-order contributions are 
reduced to products of radial matrix elements, which can be evaluated numerically. 
Sectional reports the results of such calculations for the noble gas atoms, which confirm 
that the correlation corrections to the annihilation amplitude have a marked effect on 
the gamma spectra. We also demonstrate the importance of extrapolation over the 
maximal orbital angular momentum of the intermediate electron and positron states. 
In the Appendix, expressions for the annihilation amplitudes involving many-particle 
wavefunctions in coordinate space are given. 

2. Basic Theory 

2.1. Annihilation operator 

Due to the conservation of momentum, annihilation of an electron-positron pair 
must lead to the emission of at least two photons§. In fact, according to quantum 
electrodynamics (QED), annihilation into two photons is possible only if the total spin 

§ When this process occurs in an external field, e.g., when the electron is bound in an atom, annihilation 
into a single quantum is also possible. However, its probability is small, see, e.g., Johnson et al 1964. 
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S of the pair is zero (Berestetskii et al 1982). For 5* = 1 the annihilation results in the 
emission of three photons||. 

In QED the process of electron-positron annihilation is described by the 2nd or 
3rd-order diagram, depending on the number of photons emitted (Berestetskii et al 
1982). The photons carry away the energy of the particles, ~ 2mc^, and have large 
momenta ~ mc, where m is the electron mass, and c is the speed of light. Intermediate 
electron/positron states in the QED diagrams have similar momenta. They are much 
larger than the typical momenta of valence and inner-shell atomic electrons (with the 
exception of Is electrons in heavy atoms with Z ~ 100), or the positron momentum 
in a typical experiment^. As a result, the annihilation amplitude is independent of 
the electron and positron momenta. Hence, the effective operator for annihilation into 
gamma quanta with the total momentum P is proportional to 

X! O^ki^ka^ki+ka.P , (1) 
ki,k2 

where (^k2) is the destruction operator of the electron (positron) with momentum 
ki (k2), and the (5-function ensures momentum conservation, P = ki + k2 (Ferrell 
1956). Using electron and positron destruction operators and (^(r) in the coordinate 
representation, 

ak = / e-'^-V^(r)rfr, (2) 

6k = ^ / e-''-^(r)rfr, (3) 

where V is the normalization volume, the annihilation operator can be re-written as 
(Ferrell 1956, Chang Lee 1957) 

a(P) = j e-^P '-^(r)^(r)dr. (4) 

This equation shows that the annihilation of a nonrelativistic electron-positron pair 
occurs when the particles are at the same point. Physically, this can be explained using 
the uncertainty principle. Indeed, the spatial separation between the points where 
the photons are produced is about h/mc. This is much smaller than the Bohr radius 
o-o = h'^ /me^^ i.e. the typical atomic size, or the de Broglie wavelength of the particles 
involved (c = 137 in atomic units, where % = m = \e\ = 1). 

In equations (HJ and (jU electron and positron spin indices have been suppressed. 
The operator can be used in this form if the total spin of the electron subsystem is 
zero. In this case the annihilation rate is equal to that averaged over the positron 
spin. The rate of annihilation into photons with the total momentum P is given by the 
squared modulus of the transition amplitude of operator (j^. The total annihilation rate 
is obtained by integration over d?P/ (27r)^ and summation over the final states of the 

II For S — Q annihilation into 2, 4, . . . photons is possible, while for S = 1 the number of photons must 
be odd (3, 5, . . .). In both cases the process with the smallest number of photons dominates. 
^ In positron annihilation with matter, even if the initial positron is fast, it quickly loses energy due 
to inelastic ionizing collisions, and most of the annihilation events involve slow positrons. 
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system. The correct absolute magnitude of the 2-photon and 3-photon annihilation rates 
can be determined by comparison with the spin-averaged positronium annihilation rates 
(Berestetskii et al 1982). In the first case the squared amplitude must be multiplied by 
vrrpC, while in the second case - by [4(7r^ — 9)/3]rgac, where tq = e^/mc^ is the classical 
electron radius, and a = /he is the fine structure constant. This shows that 3-photon 
annihilation is about 10^ times slower than 2-photon annihilation, and can often be 
neglected. 



2.2. Spectrum of photons and annihilation rates 

Let \i) be the state of electrons and the positron before the annihilation, and |/) 
- the state of — 1 electrons after the annihilation. In two-photon annihilation, the 
momentum distribution of the photons, 

Wf{V) = 7,rlc\{f\da{V)\i)'\ , (5) 

determines their energy spectrum. The total energy of the two photons is 

E^i + E^2 = '2mc^ + Ei-Ef, (6) 

where Ei and Ef are the energies of the initial and final states (not including the rest 
energy of the constituent particles), and total photon momentum is p^i + p^2 = P- 

For P = the photons are emitted in the opposite direction, p^i = — and have 
equal energies, E^i = E^2 = + ^{Ei — Ef) = E^. For P 7^ the photon energy 
is Doppler shifted, e.g., for the first photon E^i = E^ + VmccosO, where V = P/2m 
is the centre-of-mass velocity of the electron-positron pair, 6 is the angle between the 
direction of the photon and V, and we assume that V c and p^i = -E^i/c ~ mc. 
Hence, the shift of the photon energy from the centre of the line, e = E.yi — E^, is 

Pc 

e = —cose, (7) 
and the photon energy spectrum is given by 

^(6) = / Wf{P)6 (e - iPccos^) (8) 

Using polar coordinates and averaging over the direction of emission of the photon, one 
obtains: 

cJ J2\e\/c ' (27r) 

On the other hand, using Cartesian coordinates and choosing the z axis along the 
direction of the photon, one obtains from (jHl): 

w{e) = \jWf{P.,Py,2e/c)'^^. (10) 

This form shows that the energy spectrum is proportional to the probability density for 
a component of the momentum P. Equation (fTUj) allows one to link w{e) to the quantity 
measured by the angular correlation of annihilation radiation (AGAR) technique. In 



-W = t/ L, WfCP)—-^. (9) 
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AGAR one measures the small angle G between the direction of one photon and the 
plane containing the other photon"*". If the direction perpendicular to the plane is x, 
then B = Px/mc. Given that the distributions of P^, Py and Pz are identical (in an 
isotropic system), one obtains the distribution of from w{t) by a simple change of 
variable, = 2e/mc^. 

When a low-energy positron annihilates with a bound electron whose energy is En-, 
the centre of the photon spectrum is shifted by en/2 relative to mc^ . The width of 
the two-photon momentum distribution is determined by the typical momenta of the 
bound electron, P ~ (2m|£„|)^/^. The corresponding Doppler width of the annihilation 
spectrum, e ~ Pc ~ (|£:„|mc^)^/^, is much greater than its shift. Hence, one can regard 
the line as centred on = mc^ = 511 keV, even for the annihilation on inner-shell 
electrons. 

The total annihilation rate in the state \i) is obtained by integration over the 
momentum P and summation over the final states, 

a = E/w^/(p)|^ (11) 

d^p 

Using closure in the subspace of (A^ — l)-electron states, J2f I/) (/I = substituting (jH) 
and integrating over P, we have: 

A = vrrpC / (i|n_(r)n_|_(r)|i)dr, (13) 



/ 

vrr, 



ocE J (^|Oi(P)|/)(/|0,(P)|z)^. (12) 



where n_(r) = 'ijj^{r)'!^{r) and n+{r) = 0'^{r)0{r) are the electron and positron density 
operators. The annihilation rate is thus given by the expectation value of the electron 
density at the positron integrated over the positron coordinates. 

Equation ()13|1 gives the two-photon annihilation rate for a bound positron-atom 
state \i). The annihilation rate for a positive-energy positron moving through a gas of 
atoms is given by a similar formula (Eraser 1968), 

A = TcrlcnZes, (14) 

where n is the number density of the gas, and 

Zes = J {i\n^{r)n+{r)\i)dr. (15) 

Here the state \i) describes a collision between the positron with momentum k and the 
atom (usually in the ground state). It is normalised to the positron plane wave e*'*^'" at 
large positron-atom separations, i.e., to one positron per unit volume. The dimensionless 
parameter is interpreted as the effective number of target electrons which contribute 
to the annihilation. 

+ This is so-called ID-ACAR. A more advanced technique, 2D-ACAR, involves measuring the angle 
between the directions of the photons, i.e. two projections of the momentum. The corresponding 
distribution is proportional to / Wf{Px, Py, Pz)dPz, and is useful for studying the electron density and 
momentum distribution in anisotropic systems, e.g., solids (Puska and Nieminen 1994). 
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If one neglects the interaction between the positron and the atom and writes the 
initial state as \i) = VVbl^\0), where |0) is the ground state of the A^-electron atom, 
equation (fTHjl gives Zes = N. In reality, the values are Zes can be very different from 
the actual number of electrons (see, e.g., Surko et al 2005). Repulsion from the nucleus 
prevents the positron from penetrating into the atom, suppressing the probability of 
its annihilation with the inner electrons. On the other hand, the long-range positron- 
atom attraction increases the positron density near the atom, making Z^^ larger. The 
Coulomb attraction within the annihilating electron-positron pair increases Zgff further 
(Dzuba et al 1996, Gribakin and Ludlow 2004). 

Note that the rate (fT^ is related to the positron-atom annihilation cross section 
cTq by A = cTaUv, where v is the incident positron velocity, so that (Jq = vrrgZeffc/f. On 
the other hand, equation ()15|) shows that Z^w has the form of a transition amplitude 
between two identical states \i). Hence, a perturbation series expansion for this quantity 
can be developed (Dzuba et al 1993). 

3. Many-body theory 

3.1. Bases and building blocks 

When studying positron-atom interactions by many-body theory methods, it is 
convenient to use the bases of electron and positron states of the Hartree-Fock 
Hamiltonian of the ground-state atom. Denoting the corresponding electron (positron) 
wavefunctions and destruction operators by ip^ {ip^) and ipu)^ respectively, we have: 

Substituting into (j3)) we obtain the effective annihilation operator in the form, 

a(p) = EE(PI^I/^^)«A. (17) 

where 

(P|(5|/iz/) = 1 e-*P->^(r)y.,(r)dr (18) 

= I e-^P-('-^+'-^)/2^(ri - r2)^^(ri)<^,(r2)dridr2. 

The second form explains the presence of 6 in our notation for the amplitude. Note that 
the sum over fi in equation (fTTj) contains two distinct contributions, = J2ii<f + J2fj.>F- 
The first sum includes the electron states occupied in the atomic ground state, i.e., those 
at or below the Fermi level F ("holes"). The corresponding terms in the operator (fTTjl 
describe positron annihilation which leads to creation of a hole, as shown by diagram (a) 
in figure^ The second sum is over the electron states above the Fermi level ( "particles" ) , 
and describes positron annihilation with an excited electron, diagram (b) in figure 

Substituting (fT^ and their analogues for ip\r) and <^^(r) into n_(r) and n^{r), 
one obtains the electron-positron contact density operator from equation ()13|1 or (|T5|l . 
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Figure 1. Annihilation amplitude for a positron in state e and electron in state n 
occupied in the atomic ground state (a), or an excited electron in state fi (b). The 
double dashed line corresponds to the annihilation gamma quanta with momentum P. 



in the form, 

J n_(r)n+(r)dr = ^ ^(z/V|5|/iz/)aJ,,a]^,a^a,„ (19) 

where /x, fi' and z/, u' are the electron and positron states, respectively. The amplitude 

{iy'fi'\6\H = J ^:,{r2)rA^i)6{rr - r2)^^(ri)<^,(r2)dridr2. (20) 

corresponds diagrammatically to a point-like S vertex with two positron and two electron 
(particle or hole) lines. It is related to the annihilation amplitude (|T8|) by 

{u'^,'\6\^,u) = J (z.V|5|P)^(P|5|/.z.). (21) 

Graphically, this equation is equivalent to connecting the vertices (a) or (b) from figure 
with their complex conjugates (mirror images) and "pulling them together" by the 
double-dashed lines, to form a 6 vertex with four external lines. This basic relation 
allows one to relate products of terms in the diagrammatic expansion of the annihilation 
amplitude {f\Oa{P)\i) to the contributions to the total annihilation rate, along the lines 
of equations (fT^ and (fT^ (see below). 



3.2. Annihilation amplitudes 

Consider a process of two-photon positron annihilation with a ground-state atom, 

e+ + A — ^ + 27. (22) 

The simplest final state / of the positive ion A'^ is that of an atom with a hole in 
an electron orbital n. In the lowest order the amplitude of this process is shown by 
diagram (a), figured The exact amplitude {f\Oa{P)\i) can be represented by a many- 
body theory expansion in powers of the residual electron-positron and electron-electron 
Coulomb interactions. Figure El shows the corresponding 0th and Ist-order diagrams, 
together with the main types of 2nd-order diagrams. 

Analytically, the first diagram in figure El is equal to {P\6\ne), cf. equation 
(^H))- Expressions for the higher-order diagrams are constructed using standard 
atomic many-body theory diagrammatic rules (see, e.g., Amusia and Cherepkov 1975), 
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A3 A4 A5 



Figure 2. Many-body theory expansion of the positron-atom annihilation amplitude 
leading to the final state with a hole n. The line starting with e corresponds to the 
positron; v and fi are the intermediate positron and electron states, respectively. Wavy 
lines correspond to the electron-positron or electron-electron Coulomb interactions. 



modified to account for tlie attractive electron-positron Coulomb interaction. Tliese 
expressions contain products of Coulomb and (5-function matrix elements divided by 
energy denominators. Each energy denominator corresponds to an intermediate state 
separating two interactions in the diagram. It is given by -E — I] ^part + Z) ^hoie + ^0, where 
the sums are over all particles* and all holes in the intermediate state, E is the energy 
entering the diagram, and the infinitesimal imaginary iO defines the way of bypassing 
the pole. Summation over all intermediate positron and/or electron states is carried 
out. The overall sign factor for a diagram is given by (^—l^^h+Ni+Nep ^ where N^, Ni 
and Nep are the numbers of internal hole lines, electron-hole loops and electron-positron 
interactions, respectively. 

For example, the Ist-order diagram Al in figure |21 corresponds to 

^ {P\6\H{i^f^\V\ne) ^ ^23) 



where 



{uii\V\ne) = I ^1{y2)%{vi)- ^ -^„(ri)^,(r2)dridr2, (24) 

is the Coulomb matrix element. Diagram Al represent a correction to the annihilation 
vertex (diagram AO) due to the electron-positron interaction. Note that for low positron 
energies the denominator in equation ()23|) is never zero (i.e., the intermediate state is 
virtual for any Sy and e^), hence iO has been dropped. 

* If the annihilation photons represented by the double-dashed line, are present in the intermediate 
state, their energy must also be included in ^ Cpart ■ Since the electron and positron energies do not 
include their rest energy m<? , the energy of the photons must be given with respect to 2mc^. 
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The 2nd-order diagrams shown in figure |21 have a clear physical interpretation. 
Thus, diagram A2 represents a higher-order correction to the annihilation vertex, similar 
to diagram Al, with an extra interaction between the annihilating electron and positron. 
Such contributions are very important because the electron-positron Coulomb attraction 
strongly enhances the amplitude of finding the two particles at the same point. Diagram 
A3 describes polarization of the atom by the positron. This polarization gives rise to 
an attractive positron-atom correlation potential which behaves as —ade'^/2r^ at large 
distances, ad being the atomic dipole polarizability. It has a large effect on the incident 
positron (see, e.g., Dzuba et al 1993, 1996). Diagram A4 is a correlation correction to 
the Hartree-Fock wavefunction of the hole. Given that the Hartree-Fock approximation 
describes atomic ground state orbitals well, one may expect such corrections to be 
relatively small. Finally, diagram A5 belongs to another type of vertex corrections 
in which annihilation takes place in the presence of a virtual hole-particle excitation. 
We will show below that this and similar diagrams can be neglected because of a 
specific cancellation between the contributions of certain pairs of diagrams in the photon 
momentum distribution. 

As explained at the end of section 12. 2| Zcs can also be presented as a many-body 
perturbation series. If the state \i) in equation (fT3j) is that of a positron incident 
on the ground state atom, Z^s is given by the diagrammatic expansion shown in 
figure El Each of the diagrams contains a 5-function vertex corresponding to the contact 
density operator (fT^ with matrix elements ((20)). The 0th order diagram ZO describes 
annihilation of the positron with an electron in a Hartree-Fock orbital occupied in the 
ground state. Its analytical expression, 

Y.MdIne) = J:[ |^„(r)n^,(r)|Mr, (25) 

n n 

contains contributions of all occupied electron states n. Higher order diagrams in figureEl 
represent correlation corrections to the annihilation vertex. Their analytical expressions 
are derived using the rules formulated above. For example, diagram Zl corresponds to 
_ ^ {en\6\fiu){ufi\V\ne) ^^g) 

cf. equation ()23|). Note that in figure El we do not show corrections to the external 
positron lines e, like that in diagram A3 of figure El These corrections can be easily 
included by replacing the positron state ipe in the static field of the atom by the positron 
Dyson orbital which accounts for the positron-atom correlation potential (Dzuba et al 
1996, Gribakin and Ludlow 2004). 

By analogy with equations ()ll |) - (fT^ and (fT3|) . the total annihilation rate is related 
to the momentum distribution of the photon pairs by 

^eff = E/ K/|Oa(P)K)P^. (27) 

This means that if the terms of the series in figure El are multiplied by their complex 
conjugates and integrated over P, the expansion for Zes from figureElmust be recovered. 
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Figure 3. Many-body theory expansion of the positron-atom annihilation rate Z^g. 
Diagrams Z1-Z8 are corrections to the 0th order annihilation vertex, diagram ZO. 
Corrections to the external positron lines (e) are not shown. 



Indeed, it is easy to see, using equation (PT|) . that the contribution of diagram AO to the 
right-hand side of equation (j2Zj) (i.e., AO times its complex conjugate AO*) is equal to 
diagram ZO, after summation over n. Similarly, the product of diagrams Al and AO*, 



{en\6\P) 



(2^ 



E 



{P\S\ixu){ufx\V\ne) 



E 



{en\6\fiiy) (u^lVlne) 



(28) 



gives Zl (after summation over n), while the product of Al and Al* becomes Z3. 

However, some contributions that appear on the right-hand side of equation (j77j) 
do not match any diagrams of the Z^s expansion. In fact, their analytical expressions 
do not even have the form that would identify them with a particular diagram. For 
example, consider a product of diagrams A5 and AO* from figure |21 

{P\5\n2i'){nin2\V\nfii){iy^i\V\nie) 



{en\S\P) Y: 



i/liinin2 



(29) 



711, 



where £27 = e + En is the energy of the two gamma quanta with respect to 2mc^, given 
by the energy conservation, E^i + E^2 = ^mc^ + e + Sn- It is easy to check that the 
expression obtained upon integration of ()29p over P, can not be drawn as a diagram 
consistent with the diagrammatic rules. Hence, it also does not correspond to any 
diagram in the expansion of Zeg. 

To solve this paradox, note that equation ()27|1 contains a sum over all final states. 
Thus, besides the process shown in figure |2l where positron annihilation leads to a 
single-hole state, one must consider processes with other types of final states. Their 
formation is a result of extra electron-positron or electron-electron correlations (which 
means that their contribution to the total annihilation spectrum should, in general, be 
smaller). The simplest of such states contains two holes and one excited electron. The 
corresponding amplitude in the lowest (first) order is shown in figure |3] 

It is easy to see how these diagrams contribute to the right-hand side of 
equation ()27|1 . Thus, Bl times Bl* (integrated over P) becomes diagram Z7, Bl times 
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Figure 4. Amplitude of positron anniliilation leading to tlie final state with two holes, 
ni and n2, and excited electron /ii. For ni ^ n2, exchange diagrams with the indices 
ni and 712 swapped and an extra factor of —1, must also be considered. 



the exchange analogue of Bl* resuhs in diagram Z8, while B2 times B2* gives diagram 
Z5. At the same time, Bl times B3*, 



does not lead to any Zgff diagram after integration over P. In fact, it does not lead 
to any valid diagrammatic expression at all. However, the matrix elements in equation 
(jnn|) match those in equation ()29|) . Using the energy conservation for the process in 
figure m £ = £27 + £fj.i — £ni — , we see that the first energy denominator in 
is equal to e„ + — £:„^ — > while the first energy denominator in is equal to 
—En — e^i +eni+£n2- Heucc, the two expressions cancel term by term, as is summed 
over the hole states n, and (jHn|) is summed over the electron and hole states ui, ui and 

Note that this cancellation of the two anomalous contributions to the total 
annihilation rate, Z^s, is exact. However, it occurs even before the integration over the 
momentum of the photons. This means that these contributions can be omitted when 
calculating the photon momentum distribution and gamma spectrum. In particular, 
diagram A5 in figure |21 can be ignored. A similar cancellation also takes place for 
other diagrams representing corrections to the annihilation vertex, in which the incident 
positron line is connected with the final-state hole by means of Coulomb interactions 
and electron-hole pairs. In all of these diagrams annihilation occurs is the presence of 
one or more virtual electron-hole pairs. 

In principle, the photon spectra corresponding to different final ionic states (e.g., 
those shown in figures El and IH) could be considered separately. In this case, the 
contributions and (IHUj) are not added together, and do not cancel each other. In 
practice, though, the total energy of the two photons or the final state of the ion are 
usually not observed. The shift of the centre of the gamma line due to the difference in 
the total photon energy is negligible (as explained in section 12. 2|) and the cancellation 
described above does take place. 



E 



{n,n2\V\nfi,){en\6\P){P\6\n2u){ufi,\V\n,e) 



(30) 



un 
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4. Evaluation of the diagrams 

Let us consider a positron with momentum k which collides with a closed-sheU atom 
in the ground state, and annihilates, creating a hole in state n. The amplitude of this 
process is shown diagrammatically in figure|2l Let us denote it by A„k(P)- Omitting the 
QED factor tttqC in equation Q, we write the photon energy spectrum in this process 
as 

1 r . . ,„.,2PdPdf^p 



Integration over e gives the contribution of hole n to the total annihilation rate, 

ZM = J w^ie)de = J l^nk(P)r^. (32) 

4.1. Zeroth order 

In the 0th approximation the amplitude yl„k(P) is given by diagram AO, figure |21 

A„k(P) = {P\S\nk) = J e-^P--V^„(r)(^k(r)dr, (33) 

where 

Mr) = -Pniir)YUn), (34) 
r 

is the Hartree-Fock wavefunction of the occupied electron state in the orbital nl. The 
positron wavefunction is normalized to a plane wave, v^k(r) ~ e*^ "", and can be written 
as a partial- wave expansion (Landau and Lifshitz 1982), 

^k(r) = E exp(z5;jP,,,(r)Pjk ■ r/kr). (35) 

Here P;^ is the Legendre polynomial, and Peh{r) is the positron radial wavefunction 
with energy e = fc^/2 and orbital angular momentum li, normalized by 

P,i,{r) ~ (7rA;)-i/2 sin(fcr - ivr/i + 5,J, (36) 

where 5/^ is the phase shift. This normalization is to a 5-function of energy in Rydberg. 

Using ()34|) and ()35|) . and expanding e"*^ "" in spherical harmonics, we integrate over 
the angular variables in equations ()33|1 and ()3H1 . and sum over the electronic magnetic 
quantum number m and spin. This gives the following expression for the spectrum of 
positron annihilation with an electron from orbital nl, 

A roc 

^ni{e) = -Y^ \A^{P)\'PdP, (37) 

Cfc ^ J2|e|/c 

where 

4'^^) = (0 Q ) 3x{Pr)Pni{r)P,iAr)dr = (P||<5.||n5), (38) 



is the reduced annihilation matrix element, in which j\{z) = ^^tt/zJx^i/2{z) is the 
spherical Bessel function, A is the angular momentum carried by the photons, and the 
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notation [/] = 2/ + 1 is used. Equations (jH7|) -(|HH |) are similar to equation (III. 6) of 
Farazdel and Cade (1977), for the angular correlation function in a positron bound 
state. 



4-2. First order 

The Ist-order order correction Al, equation (j23|l . is evaluated in a similar way, using 
(I35|l as the incident positron state e. The wavefunctions of the intermediate electron 
and positron states fi and u are written as 

Mr) = lPsMr)Y,^Mn), Mr) = lPejAr)Yi^mA^), (39) 

and the Coulomb interaction V is also expanded in the spherical harmonics. 

When the 0th and Ist-order contributions are added together, the gamma spectrum 
is again given by equation (j37j) . where the amplitude is now 



In the above expression. 



(z.^||rW||n5) = X:(-l)'^''| ^, 1^ I'^MlVyWne), (41) 

is the reduced matrix element of the Coulomb interaction within the electron-positron 
pair with angular momentum A, 

A' /i W A' / 
/ V 



W^i\\Vy\\ne) = ^MlMh]{ : : , n n (^2) 



is the usual reduced Coulomb matrix element, and 

<n. = // PeMr2)PeM''l)-^Pnl{ri)PelAr2)dr,dT2, (43) 

is the radial Coulomb integral. In the amplitude ()4()|1 the sum over /i and u 
implies summation (integration) over the orbital angular momenta and energies of the 
intermediate electron and positron states. 



5. Numerical calculations and results 



In this section we present and analyse the gamma spectra obtained by using the 0th 
and Ist-order terms in the amplitude, equations (jH7|) . (jHHj) and (jlOl), for the noble gas 
atoms. 



5.1. Details of the calculations 

The calculation starts by determining the electron wavefunctions of the ground-state 
atom in the Hartree-Fock approximation. Subsequently, sets of incident positron and 
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intermediate electron and positron states in the field of the ground-state atom are 
generated. 

In experiments, the gamma spectra of noble gases were measured with thermalized 
room-temperature positrons (/c^T = 25 meV). At such low energies only the positron 
s wave contributes significantly to the annihilation signal, higher partial waves being 
suppressed by the centrifugal barrier. The shape of the annihilation spectrum is not 
sensitive to the exact value of the positron energy (as long as it remains small), and we 
calculate the spectra for the positron momentum of k = 0.05 au, corresponding to the 
energy of e = 34 meV (~ |/cbT). 

The sets of intermediate positron and electron states with orbital angular momenta 
/ = 0-8 were calculated on a uniform mesh in momentum space, starting with ko = 0.1 
au, and increasing in steps of Ak = 0.2 au, over 59 continuous spectrum wavefunctions. 
The electron and positron wavefunctions were used to compute annihilation and 
Coulomb matrix elements, and construct the annihilation amplitudes fj40|) for a range of 
momenta P. Using those, the gamma spectra for the valence and inner subshells were 
calculated from equation (jHIj). 

It is known that single-centre expansions of the annihilation amplitudes converge 
slowly with respect to the angular momenta of the electron and positron orbitals involved 
(Mitroy and Ryzhikh 1998, Mitroy et al 2002, Gribakin and Ludlow 2002). Thus, 
Gribakin and Ludlow (2002) showed that the increments in the annihilation rate upon 
increasing the maximal angular momentum from Z — 1 to /, decrease asymptotically 
as (Z -|- |)~^. In application to the annihilation spectra, this means that the quantity 
obtained by summing over all intermediate electron and positron states with / up to 
^max converges to the value for /max = oo as, 

w(e)t'— 1 ~ w(e)[°°] - ^ 1 , (44) 

/max ~l~ 2 

where A is some constantjj. We used this relation to obtain the values of w{t) = w(e)f°°l 
by linear extrapolation of w{e)^^'^'^^^ from the last two points. 

The total spectra are found by adding the contributions of the different subshells. 
They can be compared with experimental data. Annihilation rates obtained by 
integration of the spectra over e, provide information on the relative contributions of 
various subshells to the total Zgfj, and show the importance of the correlation correction 
and extrapolation over /max- 



5.2. Results for Ar 

In this section we examine in detail the results of the calculations for Ar, using it as an 
example. The next section will present a summary of the results for all noble gas atoms. 

Figure El shows the annihilation gamma spectra for the four outer subshells in Ar. 
They were obtained from equations (jn7|) and (jlUI) by including the electron and positron 

The 2nd term on the right-hand side of 1)44(1 is only the leading term, the next one being 
oc (Zmax + 5)^^! and keeping ^ next to large /max is simply a matter of convenience. 
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intermediate states with orbital angular momenta up to /max in the Ist-order term, with 
/max = 0-8. Since the spectra are symmetric, w{—e) = w{e), only the positive energies 
are shown. The Is subshell is not included, as its contribution is small and does not 
affect the results for the total spectrum. 



0.61 ' , ' , ' , ' 1 ' 1 0.15 




Energy (keV) Energy (keV) 



Figure 5. Gamma ray spectra Wni{e) for the positron annihilation on 3p, 3s, 2p and 
2s subshells in Ar, calculated using 0th and Ist-order diagrams. The graphs show how 
the result accumulates with the increase of the maximal orbital angular momentum 
^max of the intermediate electron and positron states in the Ist-order diagram. Thick 
solid curve is the result of extrapolation to Zmax 00, equation 

In the plots the lowest (thin solid) curve labelled /max = corresponds to the 0th- 
order result (for 3p and 2p), or is very close to it (for 3s and 2s), since the contribution 
of the intermediate states with the zero angular momentum is very small. For greater 
angular momenta the contributions of successive / decrease slowly, as expected. The 
spectra obtained by extrapolation to /max — > 00 (thick sohd curves) are noticeably 
larger than those obtained with /max = 8. For the two outer subshells, extrapolation 
increases the spectra by about 20% at small e. Details of the extrapolation procedure 
are illustrated by figurelHl In agreement with equation (HH), the dependence of ^[''"'''^^(e) 
on (/max + |)"^ is close to linear at large /max- 

Figure El shows that the effect of electron-positron correlations described by the 
Ist-order correction, is strongest for the 3p and 3s subshells which have small ionization 
potentials {I^p = 15.76 eV, /3s = 29.24 eV). The electrons in the inner shells are much 
more strongly bound {l2p = 249 eV, I^g = 326 eV). This makes it harder for the positron 
to perturb their motion, and reduces the effect of correlations. Nevertheless, the 1st- 
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1 1.5 

l/(/ . +1/2) 



Figure 6. Extrapolation of the spectral density W3p{e) with respect to Imax- Different 
symbols show the spectra calculated for /„iax = 0-8 and /,„ax oo, for different 
energies: full circles, e = keV; squares, e = 1 keV; diamonds, e = 2 keV; triangles, 
e = 3 keV. The values at (^max + ^)~^ = are obtained by linear extrapolation, 
equation l|44|l . from the values for ^max = 7 and 8. 



order correction leads to a noticeable increase of the annihilation signal for them as well. 
The larger binding energies also mean that 2s and 2p electrons move at greater speeds 
and have a broader momentum distribution. As a result, their gamma spectra are wider, 
as seen in figure El Because of the repulsion from the nucleus, the positron wavefunction 
has a weaker overlap with the inner-shell electron orbitals, and their contribution to the 
positron-atom annihilation rate is much smaller than that of the valence shell. 

Tables ^ and El quantify the effect of the Ist-order correlation correction on the 
shapes of the gamma spectra and annihilation rates. Table Q gives the full widths 
at half-maximum (FWHM) of the annihilation spectra for different orbitals. It shows 
that the Ist-order correction reduces the FWHM values of the spectra. This can be 
explained by noticing that in the Ist-order diagram (Al) the positron annihilates with 
an excited electron "pulled out" of the atom. The wavefunction of such excited electron 
is less localized than the bound state wavefunction. The corresponding typical electron 
momenta are lower, making for smaller Doppler shifts and a narrower gamma spectrum. 

The results shown in tableware the breakdown of the annihilation rate Zes for the 
2s, 2p, 3s and 3p subshells. Partial Z^s were obtained by integration of the corresponding 
spectral densities if„/(e), equation (j^ . According to the table, extrapolation over /max 
beyond /max = 8 increases the Zes values by about 15% for the 3s and 3p orbitals, 
and by 2-3% for the n = 2 orbitals. The final values show that adding the correlation 
correction more than doubles the Z^s value for the n = 3 shell, compared with the 0th- 
order calculation. It also increases the contributions of the 2s and 2p subshells, by 22% 
and 28%, respectively. Given their large binding energies, this is a remarkably strong 
correlation effect. In particular, it is much greater than their contribution of the n = 2 
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Table 1. Full widths at half- maximum of the gamma spectra for Ar (in keV). 







0th ^ 


~ 1st order 


Subshell 


0th order 






2s 


5.18 


5.20 


5.15 


2p 


9.42 


9.29 


9.18 


3s 


1.86 


1.84 


1.82 


3p 


2.89 


2.75 


2.71 


Total 


2.65 


2.58 


2.55 



Table 2. Contribution of the outer and inner-shell electrons to Zos in Ar. 







0th + 


1st order 


Subshell 


0th order 






2s 


0.00237 


0.00284 


0.00289 


2p 


0.00621 


0.00777 


0.00801 


3s 


0.13585 


0.23265 


0.26610 


3p 


0.60562 


1.34068 


1.57887 


Total 


0.75005 


1.58394 


1.85587 



electrons to the positron-atom correlation potential, which can be estimated from their 
contribution to the dipole polarizability of Ar. 

The size of the correlation effect for the inner-shell electrons can be compared 
with the enhancement factor for positron annihilation on hydrogen-like positive ions 
(Novikov et al 2004). This factor is defined as the ratio of Z^s calculated with an 
accurate correlated wavefunction to the value obtained using a product of electron and 
positron densities n±(r): Z^s = J n_(r)n_|_(r)dr. The latter equation would be correct 
if the total wavefunction of the system was an uncorrelated product of the electron 
and positron wavef unctions, cf. equation ()15|1 . In the many-body theory approach, the 
enhancement factor can be defined as the ratio of Z^g calculated from the expansion in 
figure ini to the value obtained from the Oth-order diagram ZO, equation pKjl (Dzuba et 
al 1996, Gribakin and Ludlow 2004). This enhancement factor describes the effect of 
short-range electron-positron correlations (Dzuba et al 1996). For the annihilation of 
s-wave positrons on the valence electrons of noble gas atoms and hydrogen, it ranges 
from about 2.5 for Ne to 6 for Xe and H (Dzuba et al 1996, Gribakin and Ludlow 2004, 
Ludlow 2003). According to Novikov et al (2004), the enhancement factors for and 
F^"*" are 1.41 and 1.20, while the electron ionization potentials in these systems are 340 
and 1100 eV, respectively. The ionization potentials of the 2s and 2p subshells of Ar 
are similar to the first of these values, and the enhancement factors for Z^g due to the 
Ist-order correlation correction in tableware of comparable magnitude, 1.22 (2s) and 
1.28 (2p). 
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In spite of the large increase of Z^s above the Oth-order result, the final total in 
table 121 (1-86) is still much smaller than the experimental Z^s = 26.77 ± 0.09 obtained 
with room temperature positrons (Coleman et al 1975). Of course, the large size of the 
Ist-order diagram means that higher-order corrections are also important, especially for 
the valence and subvalence orbitals. The main contribution here comes from electron- 
positron ladder diagrams (e.g., A2 in figureEl). Another important effect is the distortion 
of the positron wavefunction by the positron-atom correlation potential (such as that 
shown by diagram A3). As mentioned in section 13. 2| this strong nonperturbative 
effect can be accounted for by calculating the positron wavefunction from the Dyson 
equation. All-order summation of the ladder diagram series and the Dyson equation for 
the positron have been incorporated in the calculations of (Gribakin and Ludlow 
2004, Ludlow 2003), yielding good agreement with accurate calculations for hydrogen 
and experiment for the noble gases. We are currently applying a similar approach to 
the calculation of gamma spectra (Dunlop and Gribakin 2005). 

Figure m shows our final gamma spectra for the 2s, 2p, 3s and 3p subshells. Addition 
of the partial spectra leads to a smoother, almost featureless total. The contributions 
of the inner shells are very small near the centre of the line. However, they noticeably 
push the wings up at |e| > 4 keV, and contribute to a overall non-Gaussian appearance 
of the spectrum. A similar effect of core annihilation is a common feature in AGAR 
and Doppler broadening spectra of positron annihilation in solids (see, e.g., Lynn et al 
1977, Alatalo et al 1995, 1996, Asoka-Kumar et al 1996, Eshed et al 2002). 

10° 

10" 



1^10-^ 
10-' 

"■^^ -4 -2 2 4 6 

Energy (keV) 

Figure 7. Gamma spectra for 2s, 2p, 3s and 3p subshells in Ar obtained in the 0th 
+ Ist-order calculation: , 2s; , 2p; — • — , 3s; — • • — , 3p; , total. 

Let us now compare the total gamma spectrum with experiment. In Iwata 1997b, 
parameters of a large number of gamma spectra for atomic and molecular species were 
reported. The shapes of the spectra were determined by fitting the observed spectra 
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with a two-Gaussian function, 



q{E) = exp 



E — Eq 
aAEi 



+ D exp 



E — Ec 
aAEn 



(45) 



convolved with the detector response function (see also Iwata 1997c). In equation (j45p . 
a = (41n2)~^/^, and the fitting parameters are the line centre Eq, the FWHM of the 
two Gaussians, AEi and AE2, and the relative weight of the 2nd Gaussian D. The 
line-shape parameters for the noble gas atoms are given in table El 



Table 3. 7-ray line-shape parameters from two-Gaussian fits (fwata et al f997b,c). 



Atom 


AEi (keV) 


/^E2 (keV) 


D 


He 


2.15 


3.90 


0.177 


Ne 


3.14 


6.12 


0.060 


Ar 


2.25 


7.27 


0.010 


Kr 


2.02 


6.86 


0.016 


Xe 


1.80 


5.03 


0.033 



In figure IHl we compare the experimental spectrum for Ar in the form ()45p. with the 
total spectra obtained in the 0th and Oth-|-lst-order approximations. For the purpose 
of comparison, all spectra are normalized to unity at zero energy shift e. As discussed 
above, inclusion of the Ist-order correction makes the spectrum narrower. It is also 
clear that the enhancement due to this contribution is stronger near the centre of the 
line and weaker in the wings. As a result, the wings of the Oth-|-lst-order spectrum in 
figure IHl appear to be suppressed. The experimental data clearly favour the curve which 
includes the correlation correction. 




Energy (keV) 



Figure 8. Comparison of the total spectrum for Ar with experimental data. All 

spectra have been normalised to unity at e = 0. Theory: , 0th order; , 

Oth-f-lst order. Experiment: • , Iwata et al 1997b, c. 
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It must be said though, that the overall shape of the gamma spectrum and, in 
particular, its FWHM value, appear to be very robust. They are much less sensitive to 
the correlation correction than the absolute magnitude of the annihilation rate (as seen 
earlier in tables ^ and |21). The reason for such robustness of the shapes of the gamma 
spectra is that for slow positrons the momenta P of the annihilating electron-positron 
pair are determined mainly by the momentum distribution of the bound electron. Such 
distribution is described well even at the Hartree-Fock level. This explains why Iwata et 
al (1997a) were able to obtain good fits of the experimental data with crude Oth-order 
spectra, by adjusting the relative amounts of the outer and inner-shell annihilation. 

5.3. Summary of results for all noble gas atoms 

Table H] shows the FWHM values of the gamma spectra calculated using the 0th and 
Oth+lst-order approximations, as described in sect ions l5 . 1 1 and For all atoms (except 
He) the spectrum contains the contributions of the two outermost shells, e.g., for Xe, 
annihilation on 4s, 4p, 4d, 5s and 5p electrons has been considered. For heavier noble 
gas atoms the contribution of inner-shell annihilation becomes more important (Iwata 
1997a). However, it never exceeds few per cent, and has only a small effect on the 
FWHM values shown in table |3 We see that for all atoms except He, the addition of 
the correlation correction improves the agreement with experiment. 

Table 4. FWHM of the gamma spectra for the noble gas atoms (in keV). 



Theory 



Atom 


0th order 


Oth+lst 


Experiment^ 


He 


2.53 


2.35 


2.50 ±0.03 


Ne 


3.82 


3.63 


3.36 ±0.02 


Ar 


2.65 


2.55 


2.30 ±0.02 


Kr 


2.38 


2.30 


2.09 ±0.02 


Xe 


2.06 


1.99 


1.92 ±0.02 



^ Values obtained from Gaussian fits to the data (Van Reeth et al 1996, Iwata et al 
1997b). 

The changes in Zes (table Ej) due to the Ist-order correction to the annihilation 
amplitude, are more pronounced. The corresponding enhancement of the annihilation 
rate ranges from a factor of 1.9 in Ne to 2.9 in Xe. It is larger for heavier noble gas 
atoms, where the electrons are less strongly bound, making the correlation effect greater. 
Nevertheless, as in Ar, the final Zcs values fall far short of experiment. Higher-order 
corrections to the annihilation vertex ultimately increase it above the Oth-order result 
by factors of 2.5 to 6, with a further increase coming from the distortion of the positron 
wavefunction by the positron-atom correlation potential (Dzuba et al 1996, Ludlow 
2003). The latter effect becomes progressively stronger in Ar, Kr and Xe, due to the 
existence of virtual positron states with energies closer to zero (Dzuba et al 1996). 
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Table 5. Z^b values for the noble gas atoms. 



Theory 



Atom 


0th 


Oth+lst 


Experiment 


He 


0.688 


1.395 


3.94 ±0.02^ 


Ne 


0.975 


1.850 


5.99 ±0.08^ 


Ar 


0.750 


1.856 


26. 77 ±0.09^ 


Kr 


0.703 


1.870 


64.6±0.08^ 65.7 ±0.3*^ 


Xe 


0.645 


1.887 


400-450'=, 401 ± 20^^ 



^ Measured in the gas (Coleman et al 1975). 

^ Measured in the gas, with small amounts of He or H2 added to Xe to improve positron 
thcrmalization (Wright et al 1985) 
Measured in the positron trap (Murphy and Surko 1990). 

Finally, in figure IHl the gamma spectra calculated in the 0th and Oth+lst-order 
approximations are compared with the two-Gaussian fits of the experimental data (Iwata 
1997b, c). In all cases the inclusion of the vertex correction leads to better agreement 
with experiment. In particular, the relative contribution of the "shoulders" associated 
with the inner-shell contributions in Kr and Xe, is described more accurately. As in Ar, 
the enhancement of the annihilation rates due to correlations is greater for the outer up 
and ns electrons than for the inner {n — l)d, (n — l)p and {n — l)s subshells. When the 
total spectra are normahzed to unity at e = 0, this results in suppression of the wings. 

6. Summary and outlook 

In this work we have shown how a many-body perturbation series can be developed 
for the amplitude of positron-atom annihilation into photons with the total momentum 
P. We have also shown how this series converts into a diagrammatic expansion of the 
annihilation rate Zeg, upon integration over P and summation over the final states. 

Expressions for the 0th and Ist-order diagrams in the annihilation amplitude have 
been derived, in terms of reduced 5-function and Coulomb matrix elements. These 
contributions have been evaluated numerically for the valence and inner subshells of 
the noble gas atoms. We have demonstrated the importance of extrapolation over the 
angular momentum of the intermediate electron and positron states. A comparison 
of the calculated and measured gamma spectra has confirmed that inclusion of the 
correlation correction improves agreement with experiment. The correlation correction 
also leads to a sizeable enhancement of the annihilation rates, especially for the valence 
and subvalence orbit als. 

The large role played by the Ist-order diagram means that higher order 
contributions must be considered. Thus, one needs to include vertex corrections 
containing the electron-positron ladder diagrams to all orders. One also needs to 
improve the incident positron wavefunction by taking into account the positron-atom 
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Energy (keV) Energy (keV) 

Figure 9. Comparison of tiie total gamma spectra for He, Ne, Kr and Xe witti 

experiment. Tlie spectra tiave been normalized to unity at e = 0. Theory: , 0th 

order; , Oth+lst order. Experiment: • , Iwata et al 1997b, c. 

correlation potential. So far, this programme has been realised in the calculations of 
positron-atom annihilation (Gribakin and Ludlow 2004, Ludlow 2003). Work is under 
way to implement this for the gamma spectra. This will allow us to obtain accurate 
gamma spectra for positron-atom annihilation. In particular, a much more accurate 
determination of the annihilation fraction of the inner-shell electrons in the noble gas 
atoms will be possible. 

It is expected that our many-body theory work will provide further insights into the 
physical mechanisms which determine the shapes and intensities of the gamma spectra. 
We hope that this information, and possibly even some of the methods we are developing, 
will be useful for the studies of positron annihilation in more complex systems, such as 
molecules, clusters and condensed matter. 
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Appendix: Annihilation amplitude in the coordinate representation 

Suppose the initial state of the system is that of a positron with momentum k incident 
on a ground-state atom. The amplitude of finding the positron and one of the electrons 
at the same point and having a total momentum P in an annihilation event which leaves 
the ion in the final state /, is 

V A J 

X (5(r - ri)^k(ri, . . . , r^, r)dri . . . dr^ydr. (A.l) 

Here \E'k(ri, r2, . . . , r^r, r) is the initial wavefunction of positron coordinate r and A^ 
electron coordinates rj, and \l//(ri, . . . , r7v-i) is the wavefunction of the ion in the 
final state /. The initial state is normalized at large positron-atom separations 
according to \l'k(i"i, • • • , ^n, r) — '^'o(ri, • • • , riv)e*'*^''^, where $o is the target ground state 
wavefunction. Since the positron can annihilate with any target electron, the terms in 
the sum over i are antisymmetrized by the (—1)* factor. Formula ()A.1|) is analogous to 
equation (2) of Mitroy et al (2002). 

To obtain the total annihilation rate, the amplitude ()A.1|) must be squared and 
integrated over P, and a sum over all final states of the ion should be taken: 

" i:/i^/.(p)p|^ 

E5(r-r,) |^k(ri,...,r„,r)|^dri,...,dr„,dr. (A.2) 



^eff 

' / 

N 



i=l 



In the lowest approximation the initial-state wavefunction can be written as a 
product of the positron and target wave-functions, 

^k(ri, . . . , rjv, r) = <l>o(ri, . . . , r^)(^k(r) (A.3) 

where $0 is the wavefunction of the target in the Hartree-Fock approximation (Slater 
determinant) and v?k(r) is the incident positron wavefunction. Assuming that ^I// is also 
a Slater determinant in which a single-particle electron state n is missing, we obtain the 
annihilation amplitude 

A„k(P) = J e-*P-^^„(r)<^k(r)dr, (A.4) 

which is identical to the Oth-order expression 
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